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SUMMARY 


The minimum-mass structural efficiency curve is determined for sandwich- 
blade stiffened composite compression panels subjected to buckling and strength 
constraints. High structural efficiencies are attainable for this type of con- 
struction. A method of analysis is presented for the buckling of panels of this 
configuration which shows that buckling of such panels is strongly dependent on 
the through-the-thickness transverse shearing of the stiffener. Experimental 
results are presented from seven test specimens made of layered graphite/epoxy 
having stiffeners of sandwich construction with an aluminum honeycomb core, and 
there is good agreement between theory and experiment. 

INTRODUCTION 

Determination of the minimum mass of a stiffened structural panel con- 
structed of laminated filamentary materials to carry a specific compression 
load in the plane of the panel is an important design problem in aerospace 
applications. Stiffened panels having open-section stiffeners are generally 
easier to inspect and fabricate than panels having closed-section stiffeners. 
Conventional open-section stiffened panels are not as efficient, however, as 
closed-section stiffened panels (ref. 1). One way to increase the performance 
of open-section stiffened panels is by using sandwich-blade stiffeners. In 
this configuration the webs of blade stiffeners are composed of faces of high- 
strength material bonded to a low-density core such as aluminum honeycomb. By 
moving the high-strength material out from the blade midsurface the critical 
twisting and local buckling strains can be increased without appreciably 
increasing the total structural mass. Some of the advantage resulting from mov- 
ing the high-strength material out from the midsurface is offset by the decrease 
in stiffness due to one face shearing with respect to the other. This through- 
the-thickness transverse shearing depends on the shear stiffness of the core. 

Simple, closed-form column and local buckling formulas which treat various 
buckling modes independently are adequate for the design of closed-section con- 
figurations, such as the hat-stiffener configuration (ref. 2). More sophisti- 
cated buckling analyses are required, however, to predict the buckling behavior 
of open-section configurations due to their more complicated stiffener twisting 
and rolling behavior (ref. 1) . For a conventional stiffened panel acting as a 
column, inplane shearing stiffness of the webs of the stiffeners determines the 
overall transverse shearing stiffness of the column. For a stiffened panel 
with sandwich webs through-the-thickness shearing stiffness of the webs as well 
as inplane shearing stiffness of the webs must be considered. In this regard 
finite-element analyses are available which include all the effects needed. 

Such analyses, however, typically have execution times that are longer than 
desirable for the many calculations required in a structural efficiency study. 
Accordingly, to fill this need, a simplified, efficient, and reasonably accu- 
rate analysis has been derived and is presented in this paper for buckling of 
laminated sandwich-blade stiffened panels. 



In addition, the paper contains a design procedure based on the newly 
derived analysis which is used to determine the minimum mass necessary to 
carry a prescribed ratio of compressive load to panel length for graphite/epoxy 
sandwich-blade panels, and results are presented as structural efficiency plots 
and are compared with previous results. The design variables are the panel 
dimensions and ply thicknesses, and the constraints are buckling and material 
allowable strength. The validity of the analytical derivations is established 
by comparison with a finite-element structural analysis and with test results 
for seven typical panels. 


SYMBOLS 


A panel plan area 

A^ cross-sectional area of ith segment 


A 11 ' A 12 A 

A 22' A 66 f 


j 


a 44 , A55 
b 

b i 

D 11' D 12'\ 
d 22' d 6 6 J 

E 

G 

h 

I l' I 2/ I p 


J 

L 

M x ,My,M X y 

m 

N 

N X ,Ny,N X y 


inplane plate stiffnesses, see equations (A3) 

transverse shearing stiffnesses, see equations (A6) 
panel width 

width of panel segment i 

bending stiffnesses, see equations (A3) 

Young’s modulus 
shear modulus 
mesh spacing 

moments of inertia in loaded and unloaded directions and polar 
moment of inertia of beam element 

torsion constant 

panel length 

moments per unit length 

number of buckles along length 

axial stress resultant of panel 

stress resultants 
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n integer 

Qx'Qy transverse shear stress resultants 

t£ thickness of layer i 

u,v,w displacements in x,y,z directions, respectively 
W panel mass 

x axial coordinate 

Yl cross-sectional coordinate along segment i 

3 x ,By rotations of normal in web 

T warping constant 

Yxz'Yyz transverse shearing strains 

Ajj Kroneker delta, 1 if i = j or 0 if i ^ j 

e axial strain at buckling 

£ x' e y'Yxy strains 

0 angle of twist 

K x ,Ky, K X y curvatures 


MINIMUM-MASS PANEL DESIGN 
Design Procedure 

A computational procedure was developed for the design of minimum-mass 
sandwich-blade stiffened panels loaded in axial compression. Nonlinear mathe- 
matical programing techniques reported in reference 3 are used in this design 
procedure in which the cross-sectional dimensions of the panel are the design 
variables and buckling and material strength are the constraints. The buckling 
analysis is provided by the present theory and the strength analysis by maximum 
allowable strain. In design studies the panel length is fixed at 76.2 cm. In 
the present analysis it is assumed that the compression panel is simply sup- 
ported and has many stiffeners spaced along its width, so that the structural 
behavior can be modeled by a typical repeating element. The geometry and 
design variables used to represent the repeating element for the sandwich-blade 
configuration are shown in figure 1. The cross section is composed of four 
characteristic elements: (1) the skin between two adjacent attached flanges, 

(2) the skin and attached flange element, (3) the sandwich web, and (4) the 
high-stiffness cap. The sandwich web and cap constitute the stiffener. Seven 
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thicknesses and four width dimensions constitute the design variables for the 
sandwich-blade configuration with lamina ply orientation angles of 0°, 90°, 
and ±45° permitted. Orthctropic stiffnesses are assumed for all elements, and 
the additional coupling effects introduced by nonsymmetric layup patterns are 
ignored . 


Design Results 

The design procedure was used to determine minimum-mass designs required 
to carry a specified load for the sandwich-blade stiffened configurations con- 
structed of graphite/epoxy and aluminum-honeycomb material. Material proper- 
ties used in these studies are presented in table I. Thicknesses were assumed 
to be continuous variables (that is, an integer number of plies was not 
required), and neither thickness nor cross-sectional dimension design variables 
were constrained to minimum limits. To satisfy the strength requirement a max- 
imum allowable axial strain (in direction of load) of 0.008 was specified. 

The structural efficiency curve shown in figure 2 was developed by first 
generating with the design procedure seven minimum-weight designs in the load 
index (N/L) range from 0.3 MPa to 7 MPa. Arbitrary honeycomb-core densities 
ranging from 24 kg/m 3 to 160 kg/m 3 were permitted in determining the minimum 
mass of a panel required to carry a specified compression load. Transverse 
shear modulus properties for the aluminum honeycomb were taken from reference 4 
and are listed in table 1(b). It was found that for minimum mass the core den- 
sity varied from 24 kg/m 3 for lightly loaded (N/L = 2,5 MPa) to 160 kg/m 3 for 
heavily loaded (N/L = 7 MPa) panels. 

The structural efficiency curve for the solid graphite/epoxy blade config- 
uration presented in figure 2 is from reference 1 in which transverse shear 
effects were not included . The results presented in figure 2 for the hat- 
stiffened configuration are taken from reference 5 and represent slightly 
higher material stiffness properties than were used in the present study. 

Comparisons of the results shown in figure 2 show the graphite/epoxy 
sandwdch-blade configuration to be approximately 30 percent more efficient 
than the graphite/epoxy solid-blade configuration (that is, the sandwich-blade 
conf iguration requires 30 percent less mass than the solid-blade configuration 
to carry the same load) and 15 percent less efficient than the graph ite/epoxy 
hat-stiffened configuration. Comparing the graphite/epoxy sandwich-blade con- 
figuration results with results presented in reference 6 for aluminum-blade 
and hat-stiffened panels shows that the graphite/epoxy sandwich-blade config- 
uration provides a 60-percent mass savings over the aluminum-blade design and 
a 50-percent mass savings over the aluminum-hat design. 

A separate investigation was carried out to check the effect on the struc- 
tural efficiency of blade web transverse shearing for the solid graphite/epoxy 
blade configuration studied in reference 1. By assuming the transverse shear 
modulus only one-fourth of the shear modulus (see table I) and by assigning 0° 
oriented graphite/epoxy material properties for the core, an efficiency curve 
was obtained using the present approach. The present results obtained for the 
graphite/epoxy solid-blade configuration are only about 3 percent heavier than 
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results reported in reference 1 in which through-the-thickness transverse shear 
effects were not included. 


BUCKLING ANALYSES 
Description 

Accurate f ini te-element analyses typically have execution times that are 
longer than is desirable for the many calculations required in a structural 
efficiency study. Consequently, a rapid simplified buckling analysis at least 
an order of magnitude faster than available finite-element procedures has been 
developed. This analysis, capable of addressing the complex twisting and roll- 
ing modes of buckling behavior of sandwich-blade stiffened panels, is presented 
in appendix A. The analysis predicts panel buckling loads, considering a vari- 
ety of buckling modes that include column buckling, stiffener twisting and 
rolling, and local buckling, as well as the interaction between these modes. 

A second analysis is presented in appendix B and applies only to local skin 
buckling where the blade stiffener does not deform. The analysis in appendix B 
is a more accurate representation for local skin buckling than that given in 
appendix A. 

Buckling analyses developed in the appendixes consider the problem of sim- 
ply supported sandwich-blade stiffened wide panels loaded in axial compression. 
It is assumed in the development of these analyses that the buckling of a wide 
mult is tiff ened panel may be analyzed by studying the buckling response of a 
single repeating element of the panel cross section. The repeating element is 
similar to that shewn in figure 1 except that the skin and attached flange have 
a common center line as indicated by the sketch in appendix A. For the general 
study in appendix A the displacement functions assumed restrict the deformations 
to a mode antisymmetric about the blade stiffener. For the local buckling mode 
study in appendix B, however, the deformations may be symmetric. Results using 
the theory of appendixes A and B are compared later in this paper with other 
theory and with experiment. 

In the present analysis all through-the-thickness transverse shearing in 
the web in the unloaded direction is completely constrained but transverse 
shearing in the loaded direction is allowed. A separate analytical study not 
presented here showed that the buckling load of the panels investigated was not 
significantly affected by the assumption that the transverse shear stiffness in 
the unloaded direction is infinite. 


Analytical Comparisons 

To determine whether or not the present results are sufficiently accurate 
for the important modes of buckling and for various levels of transverse shear- 
ing the present computed approximate results are compared to other computed 
results for two typical panels. These computations are also used to indicate 
the significance of the transverse shearing effects. In a later section of 
this paper these computed results are compared with experimental results. 


5 



Typical designs for comparative study,- Two designs, which are designated 
as moderately and lightly loaded, are considered which are near minimum-mass 
proportions necessary to satisfy buckling and strength constraints for a 
76.2-cm-long simply supported panel. The cross-sectional dimensions and lam- 
inate descriptions of the panel designs considered are shown in figure 3. (For 
bracketed expressions defining laminate layups in fig. 3, see ref. 7. HC refers 
to aluminum honeycomb.) The panels are of laminated graphite/epoxy construction 
except for the web of the blade stiffener. The web is of sandwich construction 
with an aluminum honeycomb core faced with plies of ±45 graphite/epoxy. The 
high-axial-stiffness graphite/epoxy region at the outstanding cap of the blade 
stiffener (away from the skin) is assumed to act as a beam attached to the blade 
web. Although the panel configuration shown in figure 3 has eccentricities in 
the skin to permit a smooth outer surface, these eccentricities of the attached 
flange are ignored in the analysis. The panel designated as moderately loaded 
is designed to carry 1.36 MN/m when the web is constructed with a honeycomb core 
of density 160 kg/m^ and to carry 1.15 MN/m when constructed with a 48 kg/m^ 
core density. The lightly loaded panel is designed to carry 0.53 MN/m when con- 
structed with a 48 kg/m^ honeycomb-core density. The buckling strains, which 
are based on the derivations given in appendixes A and B for the two different 
designs, and from which comparisons are made between the different modes of 
buckling, are presented in figures 4 and 5 as a function of buckle length, L/m. 

Results for moder ately and lightly loaded designs.- Results for the moder- 
ately loaded design plotted in figure 4 show, for the interim region of buckle 
length (from 7 cm to 70 cm), that the buckling strain is strongly dependent on 
core density and, therefore, strongly dependent on the through-the-thickness 
transverse shear stiffness of the stiffener: for example, for a 30-cm-long 

buckle length the 48 kg/m^ density core reduces the buckling strain by 40 per- 
cent from the solution for the infinite transverse shear stiffness case. Simi- 
larly, for the same buckle length a panel with a 160 kg/m^ core reduces the 
buckling strain by 18 percent. In this interim range buckling deformations 
involve twisting and rolling of the stiffener. In the range beyond 70 cm the 
buckling strain is not as strongly dependent on cere density. In the range of 
longer buckle lengths, column effects couple with some twisting of the stiff- 
ener; then, with still longer lengths column buckling couples with conventional 
column transverse shearing (shearing through the depth of the stiffener but not 
through the stiffener thickness). In the range below 7 cm the buckling mode is 
local with only the skin and the attached flange deforming, and the present 
results are calculated using the analysis of appendix B. Of course, with the 
stiffener not deforming there is no effect of core density of the stiffener in 
this range. 

Results for the lightly loaded design plotted in figure 5 have ranges sim- 
ilar to those indicated for the moderately loaded design. Again for the interim 
region of buckle length (from 7 cm to 70 cm) the buckling strain is strongly 
dependent on core density and, therefore, on through-the-thickness transverse 
shearing of the stiffener: for example, for a 30-cm-long panel the 48 kg/m^ 

density core reduces the buckling strain by approximately 25 percent from the 
solution for the infinite transverse shear stiffness case. 

In each of the design studies the honeycomb core was oriented with the 
stiffer direction parallel to the load axis of the panel (that is, along the 
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length of the panel; see table I). Reversing this orientation was found to 
decrease the calculated buckling strain in the interim stiffener-twisting- 
mode/stif fener-rolling-mode range by approximately 15 percent- To achieve 
maximum load, therefore, it is recommended that the core be oriented so that 
its stiffest direction is along the load axis. 

Comparisons with BUCLASP2 . - When the honeycomb core is very denser trans- 
verse shearing is suppressed and the blade stiffener web may be treated as a 
conventional plate. An exact-plate-analysis BUCLASP2 is available in refer- 
ence 8 for the buckling of an assembly of conventional plates and beams which 
includes inplane shearing but does not allow through- the-thickness transverse 
shearing deformations. Comparison of present results and the exact results 
using the method of reference 8 (figs. 4 and 5) shows that the present results 
for infinite core density follow the same trends as exact results and agree to 
within a few percent throughout most of the buckle length range. In the range 
below 7 cm when local buckling occurs there is better agreement for moderately 
loaded panels (fig. 4) than for lightly loaded panels (fig. 5) because the 
inplane deformation of the stiffener that occurs when the faces are very thin 
is not modeled by the present local buckling analysis. 

Comparisons with NASTRAN^- . - In order to assess the accuracy of the present 
analysis to account for through- the-thickness transverse shear effects compari- 
sons have been made using the finite-element analysis code NASTRAN (ref. 9). 

Use of standard QUAD1 elements in NASTRAN allowed for orthotropic plate proper- 
ties and isotropic transverse shearing. Comparison of results in figure 4 shows 
that the methods of references 8 and 9 agree for infinite core density. For 
finite core density there is also very good agreement between present results 
and the results obtained using NASTRAN. 


EXPERIMENTAL RESULTS 
Specimen Description and Test Technique 

The moderately and lightly loaded sandwich-blade stiffened compression 
panel designs studied analytically in the previous section were also studied 
experimentally. A summary of design load and mass of the panels is given in 
table II, and the cross-sectional dimensions and laminate descriptions for 
these two designs are presented in figure 3. The specimens were tested in a 
two-stiffener-wide configuration as shown in figure 3 with side edges unsup- 
ported during the test. The specimen ends were potted in epoxy approximately 
2.5 cm thick, were ground flat and parallel, and were tested "flat ended" for 
a nearly clamped-end test condition. Panels were tested in lengths ranging 
from 17.9 cm to 76.2 cm. 

Specimens were constructed using 7.6-cm-wide pre impregnated tape of 
graphite in epoxy resin and aluminum honeycomb. Mechanical properties for 
these materials are presented in table I. Moderately loaded specimens were 


^NASTRAN: Registered trademark of the National Aeronautics and Space 

Administration. 
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constructed using a honeycomb density of either 48 kg/m 3 or 160 kg/m 3 . Lightly 
loaded specimens used a 48 kg/m 3 density core. The concentration of 0° ori- 
ented plies at the stiffener cap was laid up but not cured in a separate opera- 
tion and was transferred to the stiffener as a unit. Epoxy was used to fill the 
honeycomb core for approximately 0.38 cm at the junctures where the honeycomb 
meets the skin and cap. A separate adhesive layer was positioned in the web 
between the honeycomb and the graphite epoxy. The entire assembly was vacuum- 
bagged and cured in an autoclave. Adhesive and epoxy materials for the core 
were not accounted for in the calculations and resulted in panels that were 
heavier than was estimated. Representative measurements, which are presented 
in table II, indicated that the panels were 5 to 33 percent heavier than calcu- 
lated. The lightly loaded panels had more adhesive materials for the core than 
the moderately loaded panels on the basis of percent of total weight. Other 
measurements of the panels indicated very close to nominal dimensions. Experi- 
mental buckling strains were obtained from strain reversal data, and the mode 
shapes were determined from displacement-gage data and moire fringe patterns 
(ref. 10). 


Experimental-Analytical Comparison 

A comparison of experimental buckling strains and mode shapes between the 
results obtained using NASTRAN and the present theory is presented in table III. 
Graphical comparisons of the buckling strain as a function of the panel length 
are presented in figures 6 and 7. Analytical results using the analysis devel- 
oped herein are for a simply supported infinitely wide multistiffener configura- 
tion while the NASTRAN results and experimental data are for the clamped two- 
stiffener configurations shown in figure 1. 

Moderately loaded design.- Good agreement for buckling strains was achieved 
between both theories and experiments as shown in figure 6 and table III. Each 
of the 48 kg/m 3 honeycomb core density specimens buckled in a mode involving 
stiffener twisting and buckling of the skin as predicted by theory. Increasing 
the honeycomb core density from 48 kg/m 3 to 160 kg/m 3 increased the buckling 
strain by approximately 40 percent and changed the mode shape from m = 2 to 
m = 1 as predicted by theory. The buckling deformation pattern for the skin 
segments of the test panels may be interpreted using the moire fringe photo- 
graphs shown in figure 8. Back-to-back strain-gage data and deflection-gage 
data, the results of which are not presented, helped to indicate that stiffener 
twisting and rolling were occurring. In recording mode shapes the buckle of the 
skin between stiffeners at the panel ends has been ignored for long wavelength 
buckling since this condition is peculiar to the clamped-end test restraint 
(ref, 1) and the corresponding clamped-end conditions imposed on the NASTRAN 
solution. Good correlation for mode shape was obtained between NASTRAN and 
experimental results. 

Lightly loaded design .- Good agreement for buckling strains was obtained 
between the results from the present theory and the experimental results 
obtained for both of the 40-cm-long lightly loaded test panels as shown in 
figure 7 and table III. However, both panels buckled in a local skin mode 
(m = 7) rather than in a stiffener-twisting/skin-buckling mode (m = 2) as pre- 
dicted by the present theory. (See table III.) A moire fringe photograph of 
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the buckle pattern is presented below the curve in figure 7. It is not clear 
why the experimental mode shapes do not agree with theory for the lightly loaded 
design while they did agree with theory for the moderately loaded design- From 
a study of figure 4, however, the minimum strain for short wavelength buckling 
(L/m < 7 cm) is close to the minimum strain for interim wavelength buckling 
(7 cm < L/m < 70 cm) for the lightly loaded panel. Thus, the determination of 
the mode shape is sensitive to the detailed test conditions. 

The ends of the specimens were all ground flat and parallel and tested 
"flat-ended" so that they are essentially clamped. The present theory assumes 
the ends are simply supported. In the interim range of buckle length where 
twisting occurs, however, the results using the present theory agree well with 
clamped experimental and NASTRAN results (fig. 4) . Because of the good agree- 
ment between clamped test results and both simple- support and clamped theory, 
it is concluded that detailed boundary conditions beyond requiring zero deflec- 
tions are apparently not important for the buckling of sandwich-blade stiffened 
panels where twisting occurs and where through- the- thickness transverse shearing 
is important. 


CONCLUDING REMARKS 

A minimum-mass study is presented for sandwich-blade stiffened compression 
panels subject to buckling and strength constraints. This study indicates that 
high structural efficiencies are attainable for open-section stiffened panels. 
The graphite/ epoxy sandwich-blade stiffened configuration is approximately 
30 percent more efficient than a graphite/epoxy nonsandwich configuration, 

50 percent more efficient than an aluminum-hat stiffened configuration, and 
60 percent more efficient than an aluminum-blade stiffened configuration. 
Analysis suitable for use in automated design is presented for the buckling of 
sandwich-blade stiffened compression panels. The buckling strain of sandwich- 
blade stiffened panels depends strongly on the transverse shear stiffness of 
the blade web core. The results of the buckling analysis which includes trans- 
verse shear deformations agree well with finite-element results. Buckling 
experiments on clamped, one-bay, graphite/epoxy panels having two sandwich- 
blade stiffeners with aluminum honeycomb core are also presented, and the 
results of the buckling analysis also agree with these experiments. Conse- 
quently, the analysis and especially the inclusion in the analysis of the 
effects of transverse shearing of the blade web core and the assumption that 
the panel consists of an infinite number of repeating bays are validated. 

The theoretical results using simple-support boundary conditions agree with 
clamped theoretical and experimental results when twisting occurs and through- 
the-thickness transverse shearing is important. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 3, 1978 
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APPENDIX A 


BUCKLING ANALYSIS OF SANDWICH-BLADE STIFFENED COMPRESSION PANELS 

Buckling of wide, simply supported, sandwich-blade stiffened panels 
loaded in axial compression is analyzed. The analysis predicts the buckling 
load considering a variety of buckling modes including column buckling, stiff- 
ener twisting and rolling, local buckling, and interaction between the modes. 
It is assumed that the buckling of a wide multistiff ened panel may be analyzed 
by studying the repeating element of the panel cross section, which is shown 
in the following sketch: 



The model is composed of six basic segments: two skin (?) and two attached 

flange plate segments (?) located symmetrically and having a common center line, 

and a sandwich web plate segment (?) connected to a beam (?) forming a blade 

stiffener. The skin, flange, web, and beam segments are numbered 1 to 4, 
respectively, and the coordinate system is indicated in the sketch. The small 
error associated with analyzing a wide panel by studying the response of a sin- 
gle repeating element of the panel is discussed in reference 1. 

Expressions are obtained for the buckling strain e in this appendix by 
the method of virtual work in conjunction with assumed displacements. The vir- 
tual work of extension, bending, and the applied strain e at buckling for the 
various plate segments may be expressed as 
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6tt 


p b i p L 

^0 ^0 


N x^ e x + Ny 6 e y + Nxy^Yxy + M x filc x + M y 5lc y 


a 12\ 


+ M xy 6K xy - e^A n - — j(w fX 6w fX + v fX Sv, x ) 


dx dy 


(Al) 


where the N's end M's are the stress resultants and moments; the e's, y xy , 
and K *s are the strains and curvatures; and w and v are def ormations that 
occur during buckling. The axial length of the panel is L, and the width of 
each segment is bj_. The strains and curvatures are given in terms of the dis- 
placements as follows: 


e x u , x 


£y = v,y 




Yxy = u ,y + v ,x 


K y w »yy > 


K X y - -2W fX y 


The stress-strain and moment-curvature relations are: 


(A2) 


N x - A^^e x + A^2^y M x - Dii^x + ^12 K y 

N y = A 22 e y + A 12 e x M y = D 22 K y + D 12 K x / 


N xy “ A 66Yxy 


M xy = D 66 K xy J 


(A3) 


In order to account for transverse shearing of the blade stiffener, the 
following work of the shearing forces must be added to the virtual work expres- 
sion for the blade web segment 


p b 3 r L 
^0 ^0 


(Qx fi Y xz + Qy<$Y yz ) 


dx dy 


(A4) 


For the web segment the curvatures are redefined to be 


'xy 


“ ® x, x 

= 6 y,y 

= ®x,y 



(A5) 
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The additional strain-displacement and stress-strain relations for the web are 


Yxz “ w , x 

Y yz = 8 y - W,y 


2x ” A 44Y X z 
2 y = A 55 Y yz 


(A6) 


For the beam the virtual work is 




EI 1 W , XX^ W , XX + EI 2 V , XX^ V , XX + E ^® f XX^® f XX + GJ ®,X^®,X + EAu f x ( ^ u f x 


- sea|w x 6w, x + v fX 6v fX + — 0, x 60 ,x 


dx 


(A 7 ) 


where EIj is the beam bending stiffness about the 73 axis, EI2 the beam 
bending stiffness about an axis in the cross-section perpendicular to the y3 
axis, E T is the warping stiffness, GJ the torsional stiffness, EA the 
extensional stiffness, I p the polar moment of inertia, and A the cross- 
sectional area of the beam. 

The total virtual work is equal to the sum of the virtual work for each of 
the six segments. By substituting the expressions (A 2 ), (A 3 ), (A 5 ) , and (A 6) 

for stress resultants and moments and the strains and the curvatures, the vir- 
tual work equation can be expressed in terms of the displacements u, v, and 
w for plate segments in the skin and in terms of u, v, w, B x , and 3y for 
the web plate segment in the blade stiffener. By requiring that the beam be 
attached firmly to the end of the web of the blade stiffener, the u, v, w, 
and 6 deformations of the beam are given in terms of the deformations of the 
web. A characteristic value problem for the lowest applied strain e (the 
buckling strain) is obtained by considering particular expressions for the dis- 
placements that satisfy geometric boundary and continuity conditions. The dis- 
placements for each of the plate segments that have been chosen have the follow- 
ing form with the following arbitrary displacement coefficients U, U3, V, W, 
W3, w 30 ' x 3' x 30 ' and Y ' plus arbitrary parameters X and y: 

Segment Displacements 


mTTX 

u = U cos 

L 




1 


v = 0 


2b l 7T Vl\ iMTx 

w = — I V + bpW + W cos I sin 

' 71 2b ly 


(A8) 
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Segment 


Displacements 


mTTx 


u = U cos 


v = 0 


w = - (V + Wb 2 - Wy 2 ) sin 


mTTx 


(A9 ) 


mTTx 

u = (U + U 3 Y 3 ) cos 


v = V sin 


mTTx 


w = 


X Y3 ( Yy3 

W 3 sin + w 30\l " cos 

b 3 \ b 3 / 


mTTx 


sm 


/ Xy 3 

3 X = X + X 3 sin + X 30 cos 

\ b 3 


YY3\ mTTx 


cos 


( X Y3 ^y3'\ mTTx 

3t 7 = W cos + Y sin ) sin 


J 


(A10) 


The displacements of the beam are given in terms of the displacements of 
segment 3 according to 

u = u(b 3 ) - (b 4 - b 3 )v /X (b 3 )^ 


v = v(b 3 ) 

w = w(b 3 ) + (b 4 - b 3 ) 3 y (b 3 ) 

e = e y (b 3 ) 


(All) 


J 


These displacements are chosen so that they satisfy geometric continuity between 
segments, simple-support boundary conditions at x = 0,L, and have zero slope 
at yj = 0 and the corresponding point on the other side of the repeating ele- 
ment. Displacement functions are chosen in the y direction that give reason- 
able results for buckling in the local, stiffener-twisting, and column modes. 
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By assuming the expressions given for w in segments 1 and 2, accuracy in the 
local (skin-only) mode has been sacrificed for simplicity, A more accurate 
solution for local skin buckling is presented in appendix B, 

Many of the displacements assumed for the general buckling analysis pre- 
sented here are zero, constant, or linear in the y direction. The skin- 
deflection expression has a trigonometric function in the y direction to 
allow local buckling and twisting. The blade displacement expressions have two 
trigonometric functions in the y direction plus arbitrary parameters X and 
Y to allow for rolling as well as local buckling and twisting. 


Because of the way the panels are constructed, it is reasonable to expect 
that in segment 3 


- W f y ~ 0 


(A12) 


at y3 = 0,b3, As a consequence of this assumption 

X ^ 

W = — W? 
b 3 


Y 

Y = — W 30 
b 3 


(A13) 


and, therefore, in effect it is assumed that the transverse shearing strain 
^yz is zero in the web. The solution is then independent of A55, and the 
number of unknown displacement coefficients is reduced from 10 to 8, 


The stability equations that result from the method of virtual work for 
this problem can be identified from 


a ll 

a 12 

a 13 


a 22 

a 23 



a 33 “ eb 33 a 34 



a 44 


Symmetric 




0 


0 

0 

0 


U 



0 


0 

0 

0 


U3 

eb 34 


0 


0 

0 

0 


V 


a 45 

- 

efc>45 

a 46 

a 47 

a 48 


w 3 


a 55 

- 

eb 55 

a 56 

a 57 

a 58 


W30 





a 66 

a 67 

a 68 

! 

X 






a 77 

a 78 

*30 







GJ 

00 

00 

*3 
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wherein what follows A = Aj^ 


(&12/&22) 


and subscripts have been added to A, 


A^j , and D^j to indicate the segment represented: 


a ll 



2 

( 2A lll b l + 2A 112 b 2 + A 113 b 3 + EA) 


fm \ 2 / 1 \ 

a 12 = ] 1“ A 113 b 3 2 + E Ab 3j 


EA 


a 22 “ A 113 b 3 3 + Eft b 3 2 j + A 663 b 3 

j (EAb 3 ) (b 4 - b 3 ) 


ITtJT 

a 23 = “ A 663 b 3 _ 

li 


m 


a 13 = “l—/ ( b 4 “ b 3 ) 


a 33 


4 2 

2 (l~) [ Dlllbl + D 112 b2 + EI 2 + EA(t>4 - b 3 ) 2 ] + A 663 b 3 


b 33 



2 

(A 3 b 3 + 2A 3 bx + 2A 2 b 2 + EA) 


a 34 



+ 


1 

“ ° 112 b 2 




b 34 = 


2 


X /rarr 
b 3 \E 



b l b 2 



a 44 



2r 


1 

- A 44 b 3 


cos 


sm 


+ El 1 sin 2 X 



X /irrrr \ 

+ 2 — ^ — j (b4 - b 3 ) sin X cos X 
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b 44 = 


m tt 


1 _ / cos X sin X\ 

- A 3 b 3 l 1 J + EA sin-* X 


+ 2 — EA(b 4 - b 3 ) sin X cos X 
b 3 


+ %) Fz 


' mir\^ /mir\^ Ay 

a 45 = \ — j A 44 F 3 + ( — ] EI X sin A(1 - cos y) + — F 4 + — F 5 


b 3 b 3 


b 45 = 


mTT 

L 


A 3 F 3 + EA sin A(1 - cos y) + — EA (b 4 - b 3 ) (1 - cos y) cos A 

b 3 


Y 

+ “ Fg 
b 3 


mH b 3 

a 46 = — —(1 - cos A) 
L X 


— ) d 123 “ a 44 
b 3 


mTT 

a 47 = ~ f 7| 

Li 


m^T 

a 48 = — F 9| 


— ] D 12 3 - A 44 
\ b 3 


X \ 

— I d 123 - a 44 

b 3 / 


mir A 

+ — f 8~ d 223 
L b 3 


A Y 

“ “ “ T - d 663 f 10 

t>3 b 3 


/ miT / mTT 

a 55 = (— j a 44 f 11 + EI 1 (1 ~ cos 

Y / mTT \ 4 h , „ 

+ 2 — EI^— J (b 4 - b 3 ) sin Y (1 - cos y) + j F 12 

b 55 = j <j^A 3 F i;L + EA (1 - cos Y) 2 + 2 “ EA (b 4 - b 3 ) sin y (1 - cos y) 


/ \ 2 
i—\ EA 

2 I P _ 

(b 4 - b 3 ) 2 + — 

sin^ Y 

\ b 3 / 

A 

i 
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mn 

( sin y \ 

Y 

a 56 = “ ~ 

a 44 [ b 3 ~ , 

+ D 3 23 — sin y 

L 

\ Y / 

b 3 


m 

a 57 = - — a 44 f 3 - 

L 



m 

— D 12 3 f S 

L 


Ay \iniT 


ItTT Y 

a 58 = - — A 44 F 13 + — F 14 
L b 3 


a 66 



j D 113 b 3 + A 44 b 3 


a 67 = 



d 113 + A 44 


b 3 

h— (1 - cos 


A) 


a 68 



f b 3 

) d 113 + a 44 ~ 


sin y 


a 77 



2 

d 113 + a 44 


f 7 + °663 


A 

b 3 


f 8 


a 78 



j d 113 + a 44 



F 10 


a 88 



+ 


a 44 f 15 + d 663^“ 



d 663 f 10 
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and 



+ 



j D 221 b l + 8 




d 663 f 8 



+ 



<b 4 - b 3 ) 2 + Erf- 
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b 3^1 - cos (X - T) 
2 X - y 
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cos (X + y) 
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b 3 

+ — (1 - cos X) 
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(mV V X 
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\L / b 3 
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f 7 


b 3 f cos X sin X 
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b 3fl - cos (X - Y) 1 - cos (X + y) 


X - Y 


X + y 


t>3f 1 - cos (X - Y) 1 - cos (X + Y) 

Ft Q + 

2 X-y X + Y 


b 3 / cos Y sin Y \ b 3 

F 11 = — 1 + J + b 3 - 2 — sin Y 


'mTT \ 2 /Y f 

~ ) d 663 f 16 + f ~ j d 223 f 15 


/ mTT /ihtt\^ 

- ,b4 ■ b3, ‘ + Er r 


+ GJf — I sin z Y 


b 3 / cos Y sin Y \ b 3 
F 33 = - — 1 + + — sin Y 

2 V Y / Y 


( d 123 f 15 + D 663 F 16> 


F n c = — 1 + 


3/ cos Y sin Y 


F lfi = — 1 - 


3 [ cos Y sin Y 


/ mfr \ X Y 

\lTj d 663 f 10 - d 223 f 9 

/mTT \ / nflT \ /irm\ 

+ E^l — j (b 4 - b 3 )^ + ETf — j + Gj|— ) sin Y cos X 
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/miT \ 2 


(r) 

(fc>4 - b 3 ) 2 + — 

A 


sin Y cos 


X 


The stability equations (A14) are a homogeneous set. To obtain a nontriv 
ial solution for the buckling strain the determinant of the matrix must vanish 
The buckling strain is the lowest value of the applied strain e that satis- 
fies this equation for various values of X, y, and m. As in reference 1, 
suggested trial values are X = 0,17r, 0,371, 0.5 tt, 0.7tt, and y = X + O.Stt, 
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LOCAL BUCKLING OF SANDWICH- BLADE STIFFENED COMPRESSION PANELS 

An analysis is developed for the skin-only local buckling of wide, simply 
supported, sandwich-blade stiffened panels loaded in axial compression* The 
method of analysis used is the method of virtual work in conjunction with cen- 
tral finite differences. For this mode of buckling the blade stiffener does 
not deform but serves to clamp the skin at the point of attachment. A sketch 
of the model studied is shown here: 


| « bi y\ b 2 


* y 

1 



For the mode considered, the slope at y = 0 is zero. For this model stretch- 
ing and bending are uncoupled and, for buckling, only bending must be consid- 
ered. The virtual work during buckling is 


Stt = 1 
o 

^»bj+b2 C ^ 

1. Jo 

M x 6K x + My6<y + 

/ *12 
K X y6K X y - eU n - ^ 


dx dy 

(Bl) 

where 

= - 

w ,xx 

"X 

M x = Dll^x + D 12 K y 





Ky = - 

w >yy 

M y = D 2 2 lc y + D 12 k x 

> 


(B2) 


K xy = 

~ 2w , xy 

M xy = D 66 K xy 






nm x 






Letting w = W (y) sin in equations (Bl) and (B2) and integrating over x 
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where A = - (a^2/a 22)* In equation (B3) the derivatives are approximated 

using central differences as indicated in the following, and the trapezoidal 
rule is used for integration: 

w n+l - 2W n + W n _ x 

(W") n = (B4 ) 

h 2 


W n+ 1 * w n 

(W* ) 1 = 

nH — h 

2 


where h is the uniform mesh spacing. Equation (B3) becomes 


(B5) 



a n <5w n 


+ 


bn^Wn+i - 2<$W n + 6w n _!) 
h 2 


A n0 + Ain' 


N=i «W n+1 

+ ^ c i 

n=0 n+— 

2 

where, away from the boundaries: 



(B6 ) 


a 


n 



°12n 


w n+l ~ 2w n + w n-l 


h 2 




w n +l 

” D 22n 


2W n + W n „i 


h2 



} (B7) 


c 1 
n+— 
2 



Wn+l " W n 


h 


J 


and is the Kroneker delta which equals zero if the subscripts are differ- 

ent and one if the subscripts are the same, and D^j and A have the sub- 


22 



APPENDIX B 


script n that refers to nodal points* The expression (B6) for the virtual 
work depends on the displacements W n specified at the nodal points within and 
on the boundaries from n = 0 to N and at the fictitious points -1 and N+l 
outside the boundaries. To apply the method of virtual work it is necessary 
that the displacements W n and the virtual displacements 6w n satisfy geo- 
metric boundary conditions. The boundary conditions specified are 

Zero slope at y = 0 or, therefore , W_], = 

Zero deflection at y = bj + b 2 or, therefore, W N = 0 

Zero slope at y = b^ + b 2 or, therefore, W^+i = 1%-1 

These conditions affect the a*s, b's, and c's, so that the following defini- 
tions must be used instead of the corresponding general expressions; 


a 0 



Wi — Wq __ 

2d 120 e ^0 w 0 

h 2 
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W X - W 0 

2d 220 

h 2 



3TDT 
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+ EAn.jWn.! 


2w N-1 + %-2 / nUT 

hN-l = ( d 22>N-1 ( D 12)N-ll~ j W N-1 
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? (B8) 




< d 66) 1 
N — 
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W N _i 

bN = 2 d 22N ~ 

h 2 
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The virtual work becomes 


Lh 

6tt = — 
2 


a 0 b! - b 0 c 1/2 \ / b 2 - 2b! + b 0 c 3 / 2 - c 1/2 \ 

) 6 W 0 + a], + ] 6w l 

.2 h 2 h / \ h 2 h 


+ , , * 


/ 


V a N-l + 


b N “ 2 b N-l + b N-2 
h 2 


C 1 “ 

c 3 \ 

N — 

N— 

2 

2 

h 


6w n _i 


(B9) 


Setting the virtual work equal to zero and, therefore, requiring that the 
coefficients of the virtual displacements vanish leads to the stability equa- 
tions (for example, with N = 8 ) : 
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(BIO) 


where 



^ d 121 + d 120 + 40 
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a 22 = D lll 



+ 4 D 121 
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h 2 
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+ (D 2 22 + 4D 221 + 2d 220)~ “ eA l — 

h 4 \ L 
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a 12 


- D 12 i + D 120 + 4D ! 

66 - 

V 2 . 



1 

2(D 2 21 + d 220)~ 
h 4 


a 78 


“ pi27 + d 126 + 40 13 

66 — 

V 2 > 



1 

2 ( D 2 27 + d 226) — 
h 4 


1 

a 13 = d 221 ~ 
h 4 


1 

a 68 = d 226 ~T 
h 4 


The D's and A's must be assigned values according to the nodal point 
represented* In the neighborhood of the jump in stiffness all the D's and 
A 1 s should be assigned values based on the segment (of width h) centered 
on the nodal points except for D55 which should be assigned a value based 
on the segment between nodal points. 

The stability equations (BIO) are a homogeneous set f and to get a nontriv 
ial solution for the buckling strain the determinant of the matrix must vanish 
The buckling strain is the lowest value of the applied strain e that satis- 
fies this equation for various values of m. Computations were based on N = 1 
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TABLE I.- MATERIAL PROPERTIES 


( a ) Graph! t e/epoxy 


Modulus in fiber direction, GPa 131 

Modulus normal to fibers, GPa 13.0 

Shear modulus, GPa 6.41 

Major Poisson's ratio 0.38 

Density, Mg/m^ . 1.52 

Thickness, mm/ply • 0.14 


(b) Aluminum honeycomb 8 


Density, 

Major transverse 

Minor transverse 

kg/m 3 

shear modulus. 

shear modulus. 


GPa 

GPa 

24 

0.117 

0.062 

32 

.159 

.083 

48 

.262 

.138 

64 

.393 

.193 

80 

.524 

.262 

96 

.641 

.324 

112 

.786 

.372 

128 

.938 

.427 

144 

00 

o 

* 

1— 1 

.469 

160 

1.21 

.483 


a Reference 4* 
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TABLE II.- DESIGN LOAD AND MASS FOR SIMPLY SUPPORTED PANELS 76.2 cm LONG 



Honeycomb- 

Theoretical 

Theoretical 

Mass/Plan 

area. 

Panel 

core 

buckling 

buckling 

kg/m 2 


density. 

load. 

strain 




kg/m 2 

MN/m 


Calculated 

Measured 

Moderately loaded 

48 

1.15 

0.0044 

6.93 

7.31 


160 

1.36 

. 0051 

7.30 

8.39 

Lightly loaded 

48 

.53 

.0046 

4.29 

5.73 


TABLE III.- BUCKLING STRAINS AND MODES 


Panel 

Honeycomb- 


Strain, 

e 

Buckles 

along panel length 

length. 

core 







cm 

density. 

Present 

NASTRAN 

Experiment 

Present 

NASTRAN 

Experiment 


kg/m 2 

theory 



theory 





Moderately loaded panel 

design 



17.9 

48 

0.0047 

0.0047 

0.0045 

1 

4 

3 

26.7 

48 

.0043 

. 0050 

.0048 

1 

a i 

a i 

30.7 

48 

.0043 

.0052 

.0042 

1 

a 2 

a i 

76.2 

48 

.0044 

.0040 

.0041 

2 

a 2 

a 2 

76.2 

160 

.0051 

-0051 

.0058 

1 

a l 

a l 



Lightly loaded panel design 



40.6 

48 

0.0046 


0-0052 

2 


7 

40.6 

48 

.0046 


-0050 

2 


7 


a Ignores buckle at panel ends resulting from clamped-end biaxial 
stress field. 
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Skin and attached flange 



0° Plies 
+45° Plies 
90° Plies 

Honeycomb 


Figure 1 Configuration options and design variables used in minimum-mass design procedure 





Load index, N/L, MPa 


Figure 2.- Compression-load structural efficiency comparison for graphite/epoxy blade, 
sandwich-blade, and hat-stiffened panels, kg/m 3 = 36.1 x 10 -6 lbm/in 3 ; MPa = 145 psi 



fii 


ffl 




(b) Lightly loaded design. 

Figure 3.- Cross-sectional dimensions of panel designs considered. 
(All dimensions are in centimeters.) 


Core density, kg/m 



Buckle length, cm 


Figure 4.- Comparison of analytical results for compressive buckling strain of 
moderately loaded sandwich-blade stiffened composite panel with panel buckle 
length. 
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Figure 7.- Compressive buckling strains for lightly loaded sandwich-blade 
stiffened composite panel. Core density, 48 kg/m 3 . 
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